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ABSTRACT 

By combining the REFLEX and NORAS cluster datasets with the NVSS radio cat- 
alogue, we obtain a sample of 145, z < 0.3, X-ray selected clusters brighter than 
3 • 10^^^ erg s~^cm~^ that show a central radio emission above 3 mJy. For virial 
masses Mvir ^ 11 clusters out of 12 (corresponding to 92% of the systems) 

are inhabited by a central radio source. This fraction decreases with higher masses as 
oc M^^^''^. If this decrease is a selection effect, it suggests that the majority of X-ray 
selected clusters host in their centre a radio source brighter than ~ 10^° W/Hz/sr. A 
division of the sample into clusters harbouring either point-like or an extended radio- 
loud AGN reveals that the steepening of the Lx — T relation for low-temperature 
clusters is strongly associated with the presence of central radio objects with extended 
jets and/or lobe structures. In the latter case, Lx oc T'^ while for point-like sources 
one recovers an approximately self-similar relation Lx oc T^'^. Monte Carlo simula- 
tions show that the steepening of the Lx — T relation is not caused by clusters being 
under-luminous in the X-ray band, but rather by overheating, most likely caused by 
the interplay between the extended radio structures and the intracluster medium. In 
the case of low-mass systems, we also find a tight correlation between radio luminosity 
and cluster temperature. The effects of the central radio source on the thermal state 
of a cluster become less important with increasing cluster mass. 

The presence of radio sources with extended structures (61, corresponding to 42% 
of the sample) is enhanced in X-ray luminous clusters with respect to 'field' radio- 
loud AGN. Furthermore, we find that the luminosity distribution of the cluster radio 
population differs from that of all radio sources, as there is a deficit of low-luminosity 
{Ln < 10^^ W/Hz/sr) objects, while the number of high-luminosity ones is boosted. 
The net effect on the radio luminosity function of radio galaxies associated to cluster 
centres is of a flattening at all luminosities La < 10^'' W/Hz/sr. 

Key words: cosmology: observations - galaxies: clusters: general - galaxies:active - 
galaxies:statistics- X-rays:galaxies:clusters - radio continuum:galaxies 



1 INTRODUCTION 

Arguably one of the most important problems in the field 
of structure formation is the question of how cooling of gas, 
which leads to star and galaxy formation, is controlled by 
various heating processes. The widely accepted, hierarchical 
picture for galaxy formation states that baryonic matter 
accretes onto dark matter halos. This picture implies 
that in the absence of any feedback process, the galaxy 
mass function follows the dark matter halo mass function. 
However, it has been known for some time that this is 
not true (see e.g. White & Frenk 1991). At the high- 



and low-mass end, there is a deficit of galaxies compared 
with dark matter halos. These discrepancies indicate 
that various heating mechanisms must be involved in the 
formation of galaxies. Galaxy clusters possess the largest 
dark matter halos found in the universe and contain gas 
at X-ray-emitting temperatures. The intracluster medium 
(ICM) of galaxy clusters has been studied intensively and 
X-ray selected clusters of galaxies provide an excellent 
laboratory to explore the complex interplay between gas 
properties and heating mechanisms. 

Galaxy clusters are known to obey fairly tight seal- 
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ing relations between global properties such as mass, 
temperature and X-ray luminosity. These relations are 
important when using galaxy clusters to estimate cosmo- 
logical parameters. It is also known that simple, self-similar 
models (Kaiser ct al. 1986) that correctly reproduce the 
scaling properties of the dark-matter distribution of galgixy 
clusters (e.g. Borgani et al. 2001; Pointecouteau et al. 2005; 
Vikhlinin et al. 2006) fail in describing the observed scaling 
relations of the hot baryonic component. This discrepancy 
increases as one moves from the regime of rich clusters of 
galaxies to that of poor clusters and groups, giving rise 
to the so-called entropy-floor problem (e.g. Ponman et al. 
1999; Ponman et al. 2003; Pratt & Arnaud 2005; Piffarctti 
et al. 2005). Some source of extra heating has to be in- 
voked to reconcile theoretical expectations and observations. 

Moreover, in the absence of non-gravitational heating, 
very dense cluster cores in which the cooling time is often 
less than the Hubble time (so-called cool core clusters) 
should cool and accrete gas at rates of hundreds and 
more solar masses per year (see Fabian et al. 2004 for a 
review). However, this model conflicts with results from 
X-ray spectroscopy that have provided evidence that there 
is no massive cooling of gas in the central regions of cool 
core clusters, (e.g. Peterson et al. 2001; 2003; Bohringer 
et al. 2002). Numerical simulations of clusters of galaxies 
that include radiative cooling (e.g. Bryan 2000; Voit & 
Brian 2001; Wu & Xuc 2002) star formation and supernova 
feedbacks (e.g. Borgani et al. 2004; 2005) are found to 
be unable to reproduce the observed properties of the 
intracluster medium. 

Currently, the most popular model that is invoked to 
explain the dearth of geis below about 1 keV (Peterson et al. 
2003; Tamura et al. 2001) in the ICM rehes on the heating 
by a central AGN. In principle, radio galaxies are energetic 
enough to halt the cooling in the centres of galaxy clusters 
and explain the high-mass truncation of the galaxy lumi- 
nosity function (see e.g. Best et al. 2006). Radio- loud ac- 
tive galactic nuclei drive strong outflows in the form of jets 
that inflate bubbles or lobes. These lobes are filled with hot 
plcisma, and can heat the cluster gas in a number of ways 
(e.g. Briiggen, Ruszkowski & Begelman 2005; Ruszkowski et 
al. 2004; Briiggen & Kaiser 2002; Churazov et al. 2001). 

Direct evidence for AGN heating is growing. Images 
taken by CHANDRA and XMM-Newton have shown that 
AGNs in the centres of groups and clusters inflate bubbles 
in the surrounding X-ray emitting gas. Observations of the 
Perseus Cluster (Fabian et al. 2003; 2006) and of the Virgo 
Cluster (Forman et al. 2005, Simionescu et al. 2006) reveal 
sound waves and weak shocks that extend out to several tens 
of Kpc. Similar results are reported from investigations of 
Hydra A (Nulsen et al. 2005), MS0735.6-I-7421 (McNamara 
et al. 2004) and Abel 478 (Sanderson et al. 2004). From a 
more statistical point of view, Croston et al. (2005) con- 
sidered a well-defined and homogeneous sample of ROSAT 
groups and report evidence for the gas properties of groups 
containing radio galaxies to differ from those associated to 
radio-quiet ones, in the sense that radio-loud groups are 
likely to be hotter at a given X-ray luminosity than radio- 
quiet groups. They attribute this to the effect of heating 
induced by radio galaxies. Dunn & Fabian (2006) analysed 



a more heterogeneous sample of 30 clusters that show clear 
signatures of radio activity either in the form of 'bubbles' 
or in the form of a central radio core and study some of the 
properties of AGN heating. 

More recently. Best et al. (2007) investigated the radio 
properties of galaxies located in the centres of 625 nearby 
groups selected from the Sloan Digital Sky Survey Data 
Release 4 (DR4; Strauss et al. 2002) and argue that AGN 
heating from the cluster central galaxy probably overcom- 
pensates the radiative cooling losses in groups of galaxies, 
therefore accounting for the observed entropy floor in these 
systems. However, Best et al. (2007) did not study the ther- 
mal properties of the ICM using X-ray data. Lastly, Jctlia ct 
al. (2006) examined the central regions of 15 galaxy groups 
with CHANDRA and find that repeated outbursts have a 
long-term cumulative effect on the entropy profiles in their 
sample. 

By following Croston et al. (2005) and Best et al. (2007), 

the present work assesses the importance of (central) ra- 
dio emission on the thermal properties of the intracluster 
medium in groups and clusters of galaxies from a statisti- 
cal point of view. Our approach provides an improvement 
with respect to previous works as it addresses the issue of 
small statistics (present in Croston et al. 2005), lack of in- 
formation on the thermal properties of the ICM (present 
in Best et al. 2007) and it also extends the investigation of 
the AGN-cluster properties up to the highest cluster masses 
(~ lO^'^M©), while both the Best et al. (2007) and Croston 
et al. (2005) analyses only go up to rich groups/poor clus- 
ters scales. Finally, we distinguish radio sources in cluster 
centres on the basis of their radio morphology. As it will be 
discussed throughout the paper, such a distinction has an 
important impact on the analysis of the general properties 
of the ICM. 

Here, we have considered a statistically significant 
sample of 145 X-ray selected clusters brighter than 3 x 
10~^^ergs"^cm"^ taken from the REFLEX and NORAS 
surveys and endowed with central radio emission above the 
3 mjy flux limit. The above limits ensure > 80% complete- 
ness both in the X-ray and radio selections. We investigate 
the thermal properties of our cluster sample and compare 
them to pre-existing results from the literature. We have also 
made use of detailed Monte Carlo simulations to study the 
behaviour of radio sources surrounded by a dense medium 
such as that of rich clusters. 

The layout of the paper is as follows: §2 presents the 
parent catalogues from which our sample was drawn, while 
§3 discusses the procedure to assign a radio source to clus- 
ter centres. §4 describes the final catalogue, with particu- 
lar attention to the issue of extended vs unresolved radio 
structures. §5 illustrates our results on the fraction of X-ray 
selected clusters which host a central radio source (§5.1), on 
the X-ray and radio luminosity relations with cluster tem- 
perature (§5.2), on the X-ray, radio and mechanical luminos- 
ity relations with the clusters velocity dispersion (§5.3) and 
on the mass deposition rate ^ (§5.4). §6 presents the re- 
sults from Monte Carlo simulations on the X-ray luminosity 
distribution of clusters inhabited by a central radio source 
and on the radio luminosity distribution of radio sources 
surrounded by a dense medium. §7 summarizes our results 
and discusses the complex interplay between radio source 
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activity and thermal properties of the clusters which host 
them. 

Throughout this work we adopt a flat cosmology with 
a matter density f2o = 0.3 and a vacuum energy density 
JIa ~ 0.7, a present-day value of the Hubble parameter in 
units of 100 km/s/Mpc, ho = 0.7. 

2 THE PARENT CATALOGUES 
2.1 NORAS and REFLEX 

The Rosat-Eso Flux Limited X-ray galaxy cluster survey 
(REFLEX; Bohringcr ot al. 2004) provides information on 
the X-ray properties, fluxes, redshifts and some identifica- 
tion details for a sample of 447, z < 0.45, galaxy clusters 
brighter than an X-ray flux of ~ 2 x 10~^'^erg s^^cm^^ (0.1 
to 2.4 KeV). It covers 4.24 sr (corresponding to ~ 34% of 
the entire sky) by observing at almost all Galactic latitudes 
|6| ^ 20° and declinations 5 ^ 2.5°. The catalogue com- 
pleteness is estimated to be > 90 per cent for fluxes brighter 
than 3 x 10^^'^crg s^^cm^'^. The number of sources fulfilling 
this flux requirement is 425. 

The Northern ROSAT All-Sky galaxy survey (NORAS; 
Bohringer et al. 2000) contains 484 clusters (including the 
supplements to the original survey catalogue) with measured 
redshifts up to 2: 0.45. This number decreases to 245 if 
one considers the flux limit for REFLEX completeness. The 
survey probes Galactic latitudes |6| ^ 20° and declinations 
5^0°, and is estimated to be ~ 82% complete with respect 
to the REFLEX survey at the same 3 x 10~^^erg s~'^cm~^ 
flux levels. 

X-ray luminosities in the [0.1-2.4] KeV band are pro- 
vided for both the REFLEX and NORAS clusters. However, 
while in the case of the REFLEX survey these are given in 
the concordance ACDM cosmology, the NORAS catalogue 
still refers to an EdS model with ho = 0.5. In order to correct 
for this problem, we converted the NORAS X-ray luminosi- 
ties from the original EdS cosmology to the one adopted in 
this work by simply writing L^°^ = ■ (xACDivi/aJEds)^, 
where x is the comoving coordinate evaluated at the speciflc 
redshift of each object. 

For the conversion from X-ray luminosities in the rele- 
vant [0.1-2.4] KeV band to cluster masses we followed the 
approach presented by Borgani et al. (2001). As a first 
step, we computed bolometric and K-corrections to the [0.1- 
2.4] KeV-observed band by using the relation provided by 
Bohringer et al. (2004) . Then, for the relation between tem- 
perature and bolometric luminosity, we considered the phe- 
nomenological expression: 

Lboi = Le (^)"(l + ^)^10*4ft-^ergs-\ (1) 

where we took Lq = 3, a — 3 (e.g. White, Jones & Forman 
1997; Wu, Xue & Fang 1999) and ^ = 0. Finally, in order to 
convert temperature into mass, we assumed the conditions 

of virialization, hydrostatic equilibrium and isothermal gas 
distribution to hold by considering the relation: 

j^ij^j [noAvi.(^)]^/^(l + ^) KeV, (2) 

(e.g. Eke et al. 1998), where we take 76% of gas to be hy- 
drogen, Mvir is the cluster virial mass expressed in solar 



masses, /3 is the ratio between the kinetic energy of dark 
matter and the gas thermal energy, and Avir(z) is the ra- 
tio between the average density p within the virial radius 
and the mean cosmic density at redshift z. In our case we 
assume /3 — 1.15, as found by the Santa Barbara clus- 
ter comparison project (Frenk et al. 1999). Once we have 
computed the virial masses, we obtain the virial radii via 
Tvir = (3Mvir/47rp)^^^. Glvcn the uncertainties associated 
with the parameters used in equations (1) and (2) and the 
errors related to the determination of Lx in both the RE- 
FLEX and NORAS surveys, we note that the values quoted 
throughout this work for both Mvir and rvir have to be con- 
sidered estimates. 

2.2 NVSS 

The NRAO VLA Sky Survey (Condon et al. 1998) is a ra- 
dio survey that has observed the entire sky north of -40 deg 
declination at 1.4 GHz. The NVSS was observed with the 
array in the D configuration (DnC for most of the southern 
sky), which provides an angular resolution of 45 arcsoc. More 
than 1.8 million sources were observed down to a flux limit 
5'i.4GHz = 2.5 mjy, and the survey is shown to be ~ 80% 
complete at fluxes brighter than 3 mJy. The rms uncertain- 
ties in right ascension and declination vary from <1 arcsec 
for relatively strong (5'i.4ghz ^ 15 m,Iy) point sources to 
~ 7 arcsec for the faintest detectable objects. 



3 MATCHING PROCEDURE 

As explained in the introduction, the scope of the present 
work is two-fold: flrst, it aims at investigating the properties 
of those radio sources that reside in the proximity of cluster 
centres. Second, it studies the impact and possible effects 
of radio sources on the surrounding ICM. Note that for our 
purposes cluster centres are assumed to coincide with the 
centres of X-ray emission, so that in most (but not all) cases 
the adopted centre is coincident with the position of the 
central cluster galaxy (hereafter CG). 

In the region where X-ray and radio surveys overlap 
{5 > —40°), we compiled a catalogue of X-ray selected clus- 
ters with enhanced central radio emission by performing 
a cross-correlation analysis between objects brighter than 
3 X 10"^^ergs"^cm"^ in the NORAS and REFLEX sam- 
ples and radio sources from the NVSS dataset with 1.4 GHz 
fluxes above 3 mJy. These flux limits were required to en- 
sure a high degree of completeness both in the primary X- 
ray selection and in the search for radio counterparts (cfr. 
Sections 2.1 and 2.2). There are 425 clusters from the RE- 
FLEX survey fulfilling the above condition on X-ray fluxes. 
This number falls to 305 when we consider only those sources 
that lie aX 5 ^ —40° , necessary for the overlapping of the ra- 
dio and X-ray surveys. The total number of clusters used for 
our analysis is then 550, where 245 come from the NORAS 
survey. 

A suitable matching radius should maximize the num- 
ber of real associations and minimize spurious matches. Here 
we have to consider three factors: i) the wide range of red- 
shifts spanned by clusters in the REFLEX and NORAS sur- 
veys (from z ~ 5 • 10~* to ~ 0.3), it) the fact that X-ray 
and radio emission in the proximity of the cluster centre can 
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be displaced from one another and in) that NORAS and 
REFLEX are flux-hmited surveys probing on average more 
luminous and more massive clusters at increasing redshifts. 
All these factors converge at indicating that we cannot use 
a fixed matching angular radius. In fact, an angular radius 
suitable for sources around z ~ 0.2 would be too large for 
the more local objects and as a result would increase (even 
by a large factor for the very local clusters) the chances of 
random associations. On the other hand, an angular radius 
suitable for local sources would most likely miss real radio- 
to-X associations at higher redshifts. 

Consequently, we have chosen to adopt a varying angu- 
lar radius in our cross-correlation analysis. Since on average 
clusters of higher mass are found to host more massive and 
therefore larger CGs (e.g von der Linden et al. 2007) - so 
that the possibility for a greater displacement between X- 
ray and radio emission originating from within the eventual 
CG is higher - we have adopted a matching radius of the 
form ^match = Q dvir , where O^ir is the angular extension of 
the cluster given by 6'vir = r^ir/x, with rvir virial radius and 
X comoving distance. 

Particular attention was paid to the determination of 
the value for the quantity Q used in ^match. We find that 
Q = 0.015 (i.e. Smatch 1.5% of the projected virial radius) 
produces a large enough sample of X-ray clusters with a 
central radio counterpart while minimizing contamination 
from interlopers. These interlopers can reside both inside 
(i.e. radio emitters other than the central source) and out- 
side (i.e. spurious matches originating from projection ef- 
fects) the clusters in our sample. 

By applying the above matching criteria to the RE- 
FLEX and NORAS samples, we find that 81 REFLEX clus- 
ters out of 305 - corresponding to 26.5% of the sample - have 
a radio counterpart that is offset from the centre of X-ray 
emission by less than 1.5% of the virial radius, while in the 
case of NORAS we have 67 radio matches out of 245 objects 
(i.e. 27.3% of the original dataset) within the same angular 
distance. It is reassuring that the percentage of matches is 
approximately the same in the two cluster surveys. These 
figures are also in agreement with the results from Best et 
al. (2007) for the occurrence of central radio counterparts in 
their SDSS cluster sample (see their Figure 2). Croston et 
al. (2005) and Dunn & Fabian (2006) find somehow higher 
values (on the order of 50%), but the selection criteria for 
these two samples are different from ours, especially since 
both groups consider systems that are on average at much 
lower redshifts than our sample. 

We note that three sources in REFLEX (namely 
RXCJ0338.4-3526, RXCJ1501.1-f 0141, RXCJ2347.7- 
2808) and two in NORAS (RXCJ1229.7+0759 and 
RXCJ1242.8-f0241) allowed for more than one radio 
counterpart within the adopted matching radius. Visual 
inspection of radio maps has shown that in all five cases 
the central radio sources had extended/multiple structures. 
The radio fiux associated with each of these objects was 
taken to be the sum of their sub-components. We will get 
back to this point in Section 4.1. 

Chances of contamination from spurious sources in the 
joint radio-X-ray catalogue as obtained above have been es- 
timated by simply vertically shifting all radio positions by 1 
degree and by re-running the matching procedure. By doing 
so, we find that in the case of both REFLEX and NORAS 
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Figure 1. Distribution of residuals between radio and X-ray posi- 
tions for objects in the REFLEX (solid line) and NORAS (dashed 
lines) catalogues. 



the probability that a radio source falls by chance within 
the search radius from the centre of X-ray emission is about 
5.6%. The adoption of a smaller matching radius sensibly 
reduces the number of radio sources found in the proximity 
of the centre of X-ray emission, while leaving the chances 
for spurious contamination almost unaltered. For instance, 
if instead of taking Q — 0.015 we consider Q — 0.01, in the 
case of REFLEX we only obtain 56 matches, while the es- 
timated percentage of random coincidences remains almost 
the same (~ 4.5%). On the other hand, increasing the al- 
lowed matching radius by even a small amount greatly in- 
creases the number of possible spurious counterparts: for 
Q = 0.02 this figure already gets as high as 13%. 

The distribution of residuals between radio and X-ray 
positions as a function of the (comoving) distance from the 
centre of the cluster is shown in Figure 1. The solid line 
indicates the case for REFLEX clusters, while the dashed 
line refers to the NORAS sample. All radio positions differ 
from the X-ray ones by less than ~ 50 Kpc. If a central 
cluster galaxy is present, this implies - as expected - that 
radio emission originates from it. 

The number of clusters brighter than 3 • 10~^^ erg 
s~^cm~^ endowed with a central radio counterpart above 
the 3 mjy level as derived from our analysis is 148. Three 
clusters, however, have multiple identifications in the NO- 
RAS and REFLEX catalogues (i.e. appear in both surveys). 
By removing these sources from the radio-NORAS sample, 
we end up with 145 clusters with a unique central radio coun- 
terpart. This is the sample that will be used throughout this 
work. 



4 DESCRIPTION OF THE SAMPLE 

Relevant properties for the radio sources found in the prox- 
imity of NORAS and REFLEX cluster centres by following 
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the procedure described in Section 3 are summarized in Ta- 
ble 2 which reports in the various columns: 

(1) Name of the parent cluster. 

(2) -(3) RA(2000) and Dec(2000) of the radio source. 

(4) Radio flux expressed in mjy (Fi.4ghz). 

(5) Log of the radio luminosity expressed in W/Hz/sr 

(logiol/fl). 

(6) Redshift of the cluster (z). 

(7) Distance between cluster centre and radio counterpart 
(Kpc). 

(8) Notes on the radio appearance, i.e. whether point-like 
(blank space) or extended/sub-structured source (denoted 
with an 'S')- 

Radio luminosities have been calculated according to 
the relation: 



100 



LlAGHz = -Fl.4GHz-D^(l + z)^ 



(3) 



where D is the angular diameter distance and a is the spec- 
tral index for radio emission [F{v) oc Since we do not 
have measurements for a, we assumed a = 0.75 which is the 
typical value found for early-type galaxies and steep spec- 
trum sources such as jets and lobes. Note that, since the 
maximum redshift probed by both cluster surveys is always 
< 0.3, estimates from equation (3) do not depend on the 
precise value of the spectral index. 

4.1 Radio sources with extended emission 

The issue of extended/sub-structured sources (point (8) in 
Table 2) deserves a digression. As mentioned in Section 3, 
five clusters are associated with radio sources that have mul- 
tiple components within the chosen search radius. A visual 
investigation of the radio images taken from the NVSS sur- 
vey for all the 145 radio objects associated with NORAS and 
REFLEX clusters shows that about 42 per cent of them dif- 
fer from a point-like structure. (Note that since the NVSS 
resolution is 45", multiple structures closer than ~ 9 Kpc at 
z = 10~^ up to ~ 250 Kpc at the highest redshifts probed by 
the REFLEX and NORAS surveys will be unresolved in the 
NVSS radio maps. We will get back to this point in Sections 
5.1 and 5.2). More precisely, this happens for 61 sources, out 
of which 33 are found in REFLEX clusters and 28 in NO- 
RAS clusters. The percentage quoted above is much higher 
than what is found for 'field' radio galaxies within the same 
redshift range and at similar radio flux levels (~ 5% cfr. 
Magliocchetti et al. 2001). We will get back to this point in 
Sections 6 and 7. 

Having labelled the radio sources with an 'S', we do 
not make any further distinction on whether the radio emis- 
sion presents extended blobs, it appears as a close or distant 
double- lobed object or features a triple/more complex struc- 
ture. For the purposes of the present work, the only division 
we are interested in is whether the radio source exhibits signs 
of extended emission or not. 

Generally, it was relatively easy to recognize different 
sub-structures that belong to the same radio source. Never- 
theless, there are a number of cases that remain uncertain 
and that have been labelled with an 'S?' in Table 2. To 
play safe, these sources have not been included in the "sub- 
structure" sample and in the following analyses they will be 
treated as point-like objects. 

Radio fluxes for extended radio sources in Table 2 have 
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Figure 2. Top panel: number of clusters per X-ray luminosity 
interval. The thick (green) lines illustrate the case for NORAS, 
while the thin (blue) ones are for the REFLEX sample. Solid lines 
show the entire Fx ^ 3 X 10~^'^erg s~^cm~^ population, while 
dashed ones only represent those objects which exhibit relevant 
(^i.4GHz ^ 3 mJy) central radio emission. Lower panel: ratio 
between number of clusters with a central radio counterpart and 
the whole Fx ^ 3 x 10~^'^erg s~^cm~^ cluster population. Open 
(green) circles are for NORAS objects, filled (blue) dots for RE- 
FLEX, while (red) squares and associated errorbars represent the 
combined REFLEX-I-NORAS dataset. 



been obtained by adding up all the fluxes of their different 
components. Radio luminosities have then been computed 
by means of equation Q for the resulting total fluxes. 



5 OBSERVATIONAL PROPERTIES OF THE 
SAMPLE 

In this section, we will show the most interesting proper- 
ties and relations of our sample. In addition to the data 
described in Section 4, some extra information such as tem- 
peratures and velocity dispersions has been gathered from 
the BAX and VIZIER databases. All the best-fit values for 
the relations presented in the following pages are summa- 
rized in Table 1, while we defer a detailed interpretation of 
the results to Sections 6 and 7. 



5.1 Fraction of radio detections 

As an interesting first exercise, we can estimate the number 
of X-ray clusters brighter than 3 erg s~^ cm^^ that host a 
central radio source above Fi.4ghz = 3 mJy. Figure [2] shows, 
as a function of X-ray luminosity, the total number of X-ray 
clusters used in this work (solid lines in the top panel of Fig- 
ure [2] the thick line represents the NORAS case while the 
thin one is for REFLEX clusters), the number of X-ray clus- 
ters with a central radio counterpart (dashed lines in the top 
panel of Figure [2]) , and the ratio between these two quan- 
tities (bottom panel of Figure (2] where squares and asso- 
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Table 1. Best-fit values for the trends presented in Section 5. All relations are of the kind logjpi; = felogjpa: + a. X-ray luminosities Lx 
are expressed in [10*'' erg/s], radio luminosities Ln in [W/Hz/sr], temperatures kT in [KeV], velocity dispersions cr in [Km/s] and mass 
deposition rates in [M0/yr]. In brackets are the number of sources used to evaluate the various fits. 
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Figure 3. Ratio between number of clusters with a central ra- 
dio counterpart brighter than 3 mjy in the NVSS maps and the 
whole Fx ^ 3 X 10~'^erg s~^cm~^ cluster population as a func- 
tion of virial cluster mass. Open (green) circles are for NORAS 
objects, filled (blue) dots for REFLEX, while (red) squares and 
associated errorbars represent the combined REFLEX+NORAS 
dataset. The dashed line indicates the best fit to the data (see 
text for details). 



Figure 4. Top panel: fraction of X-ray selected clusters with 
fluxes Fx 3 X 10~'^erg s~'cm~-^ which present central ra- 
dio emission brighter than a luminosity Lr = 10^^ ° W/Hz/sr 
as a function of virial cluster mass. Open (green) circles are 
for NORAS objects, filled (blue) dots for REFLEX, while (red) 
squares and associated errorbars represent the combined RE- 
FLEX+NORAS dataset. Bottom panel: same as above but for 
the volume-limited sample Lr ^ 10^^ W/Hz/sr, z ^ 0.14 (cfr. 
Fig.m>)- 



elated errorbars indicate the combined REFLEX+NORAS 
case) . Note that there is a clear preference for radio sources 
to appear in poor clusters. In fact, despite the relatively 
poor statistics at low luminosities, the percentage of clus- 
ters hosting a central radio- AGN brighter than 3 mJy drops 
from about [50%-100%] for Lx < 10*^ erg/s to ~ 10% at 
the highest luminosities probed by NORAS and REFLEX. 

This trend becomes even more intriguing if one consid- 
ers the fraction of clusters that present central radio emission 
above 3 mJy as a function of the cluster virial mass evaluated 
as in Section 2.1 (Figure|3]). In this case, there is a trend for 
central radio detections in clusters to decrease with cluster 
mass; we find that all but one of the 12 Mvir < IO^ ^^Mq/Io ^ 



poor clusters and groups identified in NORAS and RE- 
FLEX are associated with a central radio emitting source 
brighter than 3 mJy in the NVSS maps, while this figure 
rapidly decreases to a ~ 20 per cent as soon as one moves to 
masses Mvir > lO^'^MoV \ We can quantify this result by 
writing: fn = Nradio/Nclusters = (Afvir/M*)""", where 
M, = (3.5 ± 0.9) • 10"Mq/Io ' and a = 0.35 ± 0.08 (dashed 
line in Figure |3]) . 

Obviously, selection effects bias the results presented in 
Figures [3] and H The adopted flux limit for NORAS and 
REFLEX implies that X-ray luminous and therefore more 
massive clusters are preferentially observed at higher red- 
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shifts. On the other hand, the radio-selection function de- 
rived from having used a 3 mjy hmit for our radio sample 
(solid line in Figure 112b ) implies that while locally one is 
capable of detecting all radio sources brighter than ~ 10^^'^ 
W/Hz/sr, at the redshifts probed by massive clusters only 
radio sources brighter than ~ 10^^ W/Ifz/sr can be ob- 
served. Furthermore, probing the most luminous clusters at 
high redshifts may imply selecting relaxed clusters i.e. those 
which are cool-core ones. If radio-loud AGN are the origin 
for any feedback, then this may be an extra selection effect, 
in that cool-core clusters are more likely to host powerful 
radio sources. 

Nevertheless, the observational fact that all but one of 
the local clusters and groups host a central radio-luminous 
source suggests that - by going deep enough in radio maps - 
this might also be the case for the majority of X-ray selected 
clusters. This conclusion becomes more stringent by the fact 
that our criterion for associating radio sources with cluster 
centres necessarily misses relatively common objects such as 
triple structures where the central component is too faint in 
radio to be seen (double-lobed sources with a 'radio-quiet' 
galaxy set in their middle). 

To obtain an unbiased trend for the recurrence of radio 
sources in cluster centres, we compile a volume-limited sam- 
ple. This is presented in the bottom-panel of Figure|4]for the 
volume-limited sample z ^ 0.14, Ln ^ 10^^ W/Hz/sr (cfr. 
Figure [T2b ). chosen such that the number of sources avail- 
able for the analysis is maximised. In this case, the fraction 
of X-ray clusters with central radio emission is found to be 
around the 20 per cent level, with possibly a slight prefer- 
ence for radio sources to inhabit more massive clusters. We 
note however that this result only holds for relatively large 
clusters as the radio luminosity cut has excluded all clusters 



less massive than Mv 



0"o 



Best et al. (2007) use a sample of 625 groups and clus- 
ters of galaxies selected from the Sloan Digital Sky Survey 
and cross-correlated with the NVSS and FIRST (Faint Im- 
ages of the Radio Sky at Twenty centimeters; Becker et al. 
1995) radio samples to study the properties of radio-loud 
AGN in the brightest cluster galaxies (coinciding in the ma- 
jority of the cases with our definition of CGs, see von der 
Linden et al. 2007). In their work they show that the frac- 
tion of brightest cluster galaxies that are radio-loud AGN 
increases with (stellar) mass as oc Mg^^nj^^ up to a plateau 
level of about ~ 20 — 30 per cent reached for their highest 
mass range which is centred at around Mstoiiar ~ lO^^'^M©. 
Our sample of X-ray selected clusters allows to extend the 
Best et al. (2007) analysis to the higher-mass region probed 
by rich groups and clusters of galaxies. The trend pre- 
sented in the top panel of Figure |4] then shows the fraction 
of clusters with a central radio counterpart brighter than 
Lr 10^^'^ W/Hz/sr (limit which is comparable to the 
definition of radio loudness given in Best et al. 2005) as a 
function of cluster mass. At all cluster masses, the percent- 
age of clusters hosting a central radio-loud source is substan- 
tially constant and equal to ~ 20 per cent. Our results, both 
for what concerns the presence of a plateau in the distribu- 
tion of objects that are radio-loud AGN and for the value of 
such fraction, are then in agreement with those of Best et al. 
(2007). It is quite intriguing that the same dependence (or 
rather independence) on mass of the occurrence of central 
radio-loud sources in massive extragalactic objects is found 
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Figure 5. Fraction of radio-detected clusters that present 
extended/sub-structured radio emission as a function of virial 
cluster mass. Open (green) circles are for NORAS objects, filled 
(blue) dots for REFLEX, while (red) squares and associated er- 
rorbars represent the combined REFLEX-I-NORAS dataset. The 
dashed line indicates the best fit to the data (see text for details). 



independently of whether such a mass is that of a 'field' 
galaxy (see Best et al. 2007), of a cluster central galaxy or 
of the cluster that surrounds it. 

Another interesting point concerns the intermittency of 
the radio-loud AGN phenomenon. Best et al. (2005) argue 
that the observational finding that ~ 20-30 per cent of all 
the most massive galaxies in their sample present signa- 
tures of radio activity implies that this activity has to be 
re-triggered so often that the galaxy spends about a quarter 
of its life-time in an active state. This fuel-supply require- 
ment becomes even more demanding if one considers that 
the Best et al. (2005) selection includes radio sources that 
are brighter than those considered in this work. As already 
discussed. Figure [3] suggests that if one goes to low-enough 
radio luminosities, it is likely that the majority of clusters 
host a central radio source, just as it is found locally for our 
sample. This result would then imply an almost constant re- 
triggering of the radio source, so to allow the hosting galaxy 
to spend almost all its lifetime in an active state. 

As a final point in this section, we have considered 
those radio sources within the cluster centres that present 
extended and/or sub-structured radio emission and investi- 
gated whether there was a dependence of their recurrence 
rate on cluster mass. The results are presented in Figure [5] 
The very marginal trend for their fraction to decrease with 
Mvi. (fs = Nsub/NRADio ^ (Mvir/1.6 ■ IO'^/io^U/q)"" ') is 
most likely due to the combined effect of the resolution ca- 
pabilities of the NVSS survey with the already mentioned 
selection biases in fiux-limited surveys which preferentially 
identify massive clusters at high redshifts. This implies that, 
at variance with what happens more locally, close substruc- 
tures will not be observed as such at those redshifts, result- 
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Figure 6. Fraction of radio-detected clusters that present 
extended/sub-structured radio emission as a function of redshift. 
Open (green) circles are for NORAS objects, filled (blue) dots for 
REFLEX, while (red) squares and associated errorbars represent 
the combined REFLEX+NORAS dataset. 



ing in a net loss of sub-structured/extended sources start- 
ing from z ~ 0.05 and particularly beyond z ~ 0.1 (cfr. 
Figure Nonetheless, in the local universe and/or at low 
cluster masses there is a tendency for the majority (> 50%) 
of central radio sources to possess extended radio emission. 
Combining this result with that of Figure|3]we can then con- 
clude that, at least locally, the majority of clusters (which 
in this case are rather groups) host an extended radio source 
in their centre. We will investigate the consequence of this 
finding in the next sub-sections. 



5.2 Luminosity vs Temperature 

X-ray scaling relations of galaxy clusters, such as the X- 
ray luminosity-temperature relation, play an important role 
when using clusters of galaxies to constrain cosmological pa- 
rameters. When gravitational heating is the only source of 
heating, clusters should behave self-similarly and the X-ray 
luminosity (due to thermal bremsstrahlung) should scale as 
Lx oc (Kaiser 1986). 

Most authors (e.g. White, Jones & Forman 1997; Ar- 
naud & Evrard 1999 and more recently Popesso et al. 2005) 
report Lx ~ kT correlations with logarithmic slopes close 
to 3 in the cluster regime, though attempts to remove the 
effects of central cooling flows (Allen & Fabian 1998; Marke- 
vitch 1998) results in flatter relations. However, the above 
result is found to break at the mass/luminosity scales of 
groups of galaxies. Helsdon & Ponman (2000a, b) report a 
steeper slope of 4.9 ± 0.8 for a sample of X-ray bright loose 
groups, and 4.3±0.5 for a larger sample including both loose 
and compact groups, while Xue & Wu (2000) found a slope 
5.6 ± 1.8 from data for 38 groups taken from the literature. 
More recently, Osmond & Ponman (2004) - by analysing a 
sample of 60 galaxy groups within the GEMS project - have 



measured an Lx — kT slope that is comparable with that 
found in clusters (but see their discussion about possible 
systematics which could have affected their result). 

The observed departure from self-similarity has been 
interpreted as excess of entropy in low mass clusters due to 
some form of pre-heating generated by sources such as su- 
pernovae and/or AGN (see e.g. Tozzi & Norman 2001). Ev- 
idence for heating excess in 19 groups containing radio-loud 
AGN has recently been reported by Croston et al. (2005). 
From the theoretical point of view a number of works have 
shown that radio-loud AGN can provide a sufficiently ener- 
getic and widely distributed heating mechanism to solve the 
entropy floor problem (e.g. Briiggen & Kaiser 2001; Briiggen 
& Kaiser 2002; Bohringer et al. 2002; Fabian et al. 2003; 
Heinz et al. 2006). 

In order to investigate the effects of the presence of a 
central radio component in groups and clusters of galax- 
ies from a statistical point of view, we have gathered tem- 
peratures for 64 of the radio-detected clusters belonging 
to our sample from the BAX database. This sample was 
sub-divided according to whether the central radio source 
presented an extended (31 objects) or an unresolved (33 
objects) structure. The Lx vs kT trend for the two sub- 
populations is presented in the left-hand panel of Figure [3 
where open circles represent clusters inhabited by point-like 
radio sources, while filled squares are for clusters associated 
with central sources with extended radio emission. The dif- 
ference in behaviour between these two populations is strik- 
ing: at a given temperature, the overwhelming majority of 
clusters associated with extended radio sources have a lower 
X-ray luminosity than those that host point-like radio ob- 
jects. Or, from a different point of view, at a given lumi- 
nosity clusters associated with extended radio sources are 
systematically hotter than those that are not. 

When we make a fit to our entire sample of X-ray 
selected clusters that host a central radio source, we find a 
slope of the Lx — kT relation of fe = 2.8 ± 0.1. This slope 



becomes appreciably shallower {b = 2.35 



+0.08 \ 



if one only 



includes unresolved radio sources, and remarkably steeper 
{b = 4.0 ± 0.2) for the sub-population of extended radio ob- 
jects. Intriguingly, X-ray clusters with a point-like central 
radio source trace a distribution that is very similar to the 
Lx oc expected for self-similar systems. The presence of 
extended radio sources breaks this relation causing a much 
steeper dependence of X-ray luminosity on cluster temper- 
ature. This may point towards a more efficient mode with 
which extended sources can permeate the surrounding ICM 
and transfer heat within the cluster. As will be discussed 
in greater detail in Sections 6 and 7, the absolute effect of 
extended radio-loud AGN on the ICM decreases with clus- 
ter mass. This is clearly seen in Figure [7^. With increasing 
mass (i.e. increasing T), the difference between point-like 
and extended radio sources decreases. 

The issue of point-like vs extended radio sources has 
been investigated further by considering volume-limited 
samples of clusters. In fact, given the 45" resolution of NVSS 
(cfr Section 4.1 and Figure O, the above result on the dif- 
ferent behaviour of the Lx — kT relationship for clusters 
that host compact or extended radio sources could be bi- 
ased by the presence of radio structures that are extended 
but appear as point-like in the NVSS maps. We have then 
repeated the above analysis by splitting the samples into 
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Figure 7. Left-hand panel (a): X-ray luminosity versus cluster temperature for those X-ray selected clusters brighter than 3 erg 
cm~^ which exhibit central radio emission above 3 mjy. Open (red) circles represent clusters associated to point-like radio sources, while 
solid (green) squares are for those clusters inhabited by radio objects presenting extended structures. The dashed (green) line indicates 
the Lx — kT best fit for the sub-population of extended radio sources, while the dotted (red) one is for point-like radio objects (see Table 
1). The solid line is the best fit to the whole radio cluster sample. Right-hand panel (b): radio luminosity of the sources associated with 
cluster centres versus cluster temperature. Symbols are as before. 
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Figure 8. Left-hand panel (a): X-ray luminosity Lx versus velocity dispersion a for a subsample of 39 clusters with enhanced central 
radio emission. Open (red) circles identify point-like radio sources while filled (green) squares represent extended ones. The (green) 
dashed and (red) dotted lines respectively indicate the Lx — o" best fits to the sub-populations of extended and point-like radio sources 
(see Table 1). The solid line is the best fit to the whole radio cluster sample. Right-hand panel (b): velocity dispersion versus cluster 
temperature. Symbols are as before. 



sources with z ^ 0.05 and 0.05 < z ^ 0.2. This means 
that for the closer sub-group of clusters, all radio structures 
more extended than ^^46 Kpc will be correctly identified as 
such in the radio maps. However, in the subsample made of 
more distant objects one will identify as point-like all those 



extended radio sources that have structures smaller than 
~ 190 Kpc. In the former case, the best-fit values for the 
Lx vs kT relationship are 6 = 4.4 ± 0.2; a = -2.25 ± 0.08 
for extended sources and b = 2.29 ± 0.04; a = -1.02 ± 0.02 
for point-like structures. At higher redshifts one instead gets 
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b = 2.9^0.1; a = -1.5lo;g and b = 2.7 ± 0.3; a = -1.3 ± 0.2. 
One can then see that the results obtained for the more reli- 
able local subsets of clusters are in very good agreement with 
those derived previously in this Section (cfr Table 1) regard- 
less of the cluster redshift distribution. This is reassuring. 
Harder is the intepretation of the higher-redshift results. In 
this case clusters inhabited by point-like and sub-structured 
radio sources present a very similar dependence of the X- 
ray luminosity on the temperature. Whether this effect is 
due to the fact that the 'point-Uke' sample of sources in 
this redshift range also includes a fraction of sub-structured 
(but unresolved) objects, therefore producing a net mix of 
the two populations or whether its explanation resides in 
the more physical effect that at high redshift we are pref- 
erentially selecting higher-mass systems where the absolute 
effect of extended radio-loud AGN on the ICM is smaller is 
difficult to tell. Higher-resolution radio observations of the 
cluster centres are needed to disentangle these two effects. 

The average Lx — kT behaviour obtained by considering 
both populations of extended and point-like radio sources is 
in remarkable agreement with what is found in the litera- 
ture for the whole cluster population (e.g. White, Jones & 
Forman 1997; Arnaud & Evrard 1999; Popesso et al. 2005). 
This result, together with the fact that we expect most X- 
ray selected clusters to host a central radio object brighter 
than ~ lO'^" W/Hz/sr (see §5.1), suggests that the observed 
Lx oc scaling is due to the superposition of two separate 
populations: clusters inhabited by unresolved radio objects 
and clusters with extended radio emission. 

The above analysis underlines the fundamental impor- 
tance of extended radio-loud AGN in breaking the self- 
similarity of clusters, at least on scales kT < 3 KeV. At 
this point one can wonder how the level of radio activity 
depends on the cluster temperature. A plot of Lr vs kT is 
presented in the right-hand panel of Figure [T] Independent 
of the radio morphology, there is a much greater spread than 
in the Lx vs kT plot. However, for kT < 3 keV, i.e. in the 
regime where radio-loud AGN heating is very efficient, the 
Lr — kT relation is relatively tight, providing some evidence 
- at least in low-mass systems - for a relation between radio 
luminosity and cluster temperature, in that hotter clusters 
host more powerful radio sources. Note that this result is 
found to hold for both unresolved and sub-structured radio 
sources. 



5.3 Velocity Dispersions 

Velocity dispersion measurements for this part of our analy- 
sis have been taken from the VIZIER database and in most 
cases rely on the works of Mahdavi et al. (2000) and Mahdavi 
& Geller (2001). X-ray luminosities have been compared to 
velocity dispersions for the 39 clusters for which measure- 
ments were available in the left-hand panel of Figure [S] Best- 
fit values for the fits are provided in Table 1. The trend seen 
for the whole population of clusters inhabited by a central 
radio source {Lx oc o-', with 6 = 3.8 ± 0.4) is in general 
agreement with what is found in the literature for different 
sets of clusters and groups of galaxies (e.g. b ~ 3.6, Popesso 
et al. 2005; b ~ 4.4, Madhavi & Gheller 2001; b ~ 3.9 for the 
clusters considered in Osmond & Ponman 2004). This be- 
haviour, together with what was found in Section 5.2 for the 
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Figure 9. Ratio between mechanical and X-ray luminosity versus 
cluster velocity dispersion. Open (red) circles identify point-like 
radio sources while filled (green) squares represent extended ones. 
The (green) dashed and (red) dotted lines respectively indicate 
the imech/^x — f best fits to the sub-populations of extended 
and point-like radio sources (see Table 1); the solid line is the 
best fit to the whole radio cluster sample. 



Lx — kT relation, show that radio-identified clusters provide 
a fair sample of the entire cluster population. 

However, the similarity breaks if one considers clusters 
inhabited by unresolved and extended radio sources sepa- 
rately. For the former case, we find that b — 6.21q'3, while in 
the case of extended radio emission we find b = 5.5tJ:g. The 
difference in the observed behaviour of these two sub-groups 
with the general case (cfr. Figure |S^) may be attributed to 
the fact that in our sample extended radio sources seem 
to be systematically associated to higher velocity dispersion 
systems. Even though the available data is not good enough 
to provide more than a Icr evidence for this effect, the net 
result of combining extended and point-like radio sources is 
a flattening of the Lx — o relation to the 'standard' value 
b ~ 4. As already observed in Section 5.2, it is intriguing 
that the sum of two populations with dissimilar behaviours 
gives rise to a trend that matches previous measurements 
obtained for different cluster sets. This suggests again that 
clusters can be divided into two rather different groups: clus- 
ters that host point-like radio sources and those associated 
with extended radio emission. 

It would be interesting to investigate the Lx — o" re- 
lation in the low- mass/low- luminosity region {Lx ^ 10'*'' 
erg/s, kT < 3 KeV) where the discrepancy in the Lx — kT 
behaviour between extended and point-like sources becomes 
more pronounced. Unfortunately, the available data on ve- 
locity dispersions only probe the Lx ^ lO*^'® erg/s regime, 
so that this is currently not possible. For the same reason, 
we cannot detect any difference in the relation between tem- 
perature and velocity dispersion (right-hand panel of Figure 
[8|, a part from a marginal preference for extended radio 
sources to reside in higher velocity dispersion systems. The 
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slope of the a vs kT relation is 6 0.5 in all cases, again in 
agreement with most of the results found in the literature. 

By analysing the cavities and bubbles that are pro- 
duced in clusters and groups by the interaction between ra- 
dio sources and the surrounding hot gas, Birzan et al. (2004) 
provide a useful empirical relation to estimate the mechan- 
ical luminosity released by the central radio source into the 
ICM: 



1036W 



= (3.0 ±0.2) 



1025WHZ-1 



(4) 



Even though the 1.4 GHz synchrotron luminosity is not too 
reliable at predicting the mechanical luminosity as can be 
seen from the large errors in the exponent of Eq. U we will 
use it to calculate the energetic balance between AGN heat- 
ing and radiative losses (as in Best et al. 2007 and Birzan et 
al. 2004). This is presented in Figure |9] as a function of the 
cluster velocity dispersion, where once again solid squares 
are for extended radio sources, while open circles represent 
point-like structures. In agreement with Best et al. (2007) 
we find that, while on the smallest mass scales (cr < 400 
Km/s) radio heating balances radiative losses, such a heat- 
ing mechanism falls short by as much as a factor 100 for the 
most massive clusters in our sample. 

There is no significant difi'erence between the behaviour 
of extended and point-like sources (Lmecb/Lx oc o"~^), even 
though the data does not allow us to conclude anything 
about the group regime. 

Thus, for the smaller mass systems we have found the 
following: (i) mechanical heat and radiative losses balance 
and (ii) there is a correlation between X-ray luminosities and 
cluster temperatures and between radio luminosities and 
cluster temperatures in the low temperature regime (Sec- 
tion 5.2). This suggests the conclusion that in objects such 
as groups or small clusters of galaxies there is a strong in- 
terplay between central radio source and the surrounding 
gas. In this case, more powerful radio sources lead to hot- 
ter systems. The eventual presence of extended radio struc- 
tures such as jets and lobes is of fundamental importance 
to the thermal state of the ICM, as they can carry energy 
throughout the whole cluster and appear to be more efficient 
at heating the surrounding medium. The heating provided 
by the central radio source is also found to be sufficient to 
balance radiative losses. 

The situation is different at higher masses. In this case, 
we find that there is no correlation between the radio lu- 
minosity of the central source and the cluster temperature. 
Also, the efi'ect of the radio-loud AGN on the thermal state 
of the cluster is independent of the capabilities of such an 
object to permeate the surrounding medium (no difference 
in the Lx — kT relation between extended and point-like 
sources). AGN heating is found to be insufficient to coun- 
terbalance the radiative losses of the ICM, with a relative 
energetic importance which steeply decreases (as oc (t~^) 
with cluster mass. 

A different source of prominent heating must be invoked 
in this latter case. One such mechanism may be thermal con- 
duction (see e.g. Briiggen et al. 2003; Hoeft & Briiggen 2004; 
Roychowdury et al. 2005; Ruszkowski & Begelman 2002; Fu- 
jita & Suzuki 2005; Narayan & Medvev 2001; Voigt & Fabian 
2004). Hoeft & Briiggen (2004) have provided evidence for 
the relative importance of heating due to thermal conduc- 
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Figure 10. X-ray luminosity within the cooling volume versus 
radio luminosity Ljj of the central source. Symbols are as before. 



tion to increase with cluster mass, in agreement with our 
results. 



5.4 Cooling core clusters 

Clusters of galaxies can be separated into two classes: clus- 
ters with dense gaseous cores in which the cooling time is 
less than the Hubble time, (cool core clusters), and cluster 
with less dense gaseous cores (non-cool core clusters). 

Chen et al. (2007) have investigated the influence of cool 
cores on cluster scaling relations. Their analysis has shown 
that in cool core clusters the X-ray luminosity is enhanced 
over that of non-cool core clusters, while other parameters 
such as temperature, mass and gas mass are less affected by 
the occurrence of a cooling core. This has been explained 
by the fact that at least some of the non-cool core clus- 
ters are in dynamically young states compared with cooling 
core clusters. Alternatively, the non-cool core clusters might 
have had their cool cores heated by a radio-loud AGN. This 
is suggested by Fig. [7^ which shows the Lx — T relation for 
clusters with an extended radio source and for those without 
it. Our result indicates that the scaling relations for clus- 
ters are significantly affected by the presence of a radio-loud 
AGN at their centres. 

In the work by Birzan et al. (2004), a close correlation 
between the mechanical energy output of the jets and the 
energy loss of the central X-ray plasma by cooling (related 
to the traditionally determined cooling fiow mass deposition 
rate) has been established. This lends support to the idea 
that the feeding of the AGN is connected to the cooling rate 
and that the energy output of the AGN regulates the cooling 
core structure by feedback processes. In any model of feed- 
back, the power of the central AGN should be regulated by 
the supply of fuel from the ICM. Thus, as in a thermostat, 
large radiative losses should give rise to a strong, central 
radio source, which in turn will heat up the ICM (see Chu- 
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Figure 11. Thermal energy required to heat the clusters which 
host in their centres extended radio sources from the predicted to 
the measured temperature (cfr. Figure 7a) as a function of radio 
power (see text for details). 



razov et al. 2003). The more strongly a cluster cools at its 
centre, the more fuel can be supplied to the central AGN 
which, in turn, can drive more energetic radio lobes through 
the ICM which gets heated until the central cooling time 
goes up, the mass deposition rate goes down and the supply 
of fuel is shortened. 

Cooling flow measurements for 15 clusters belonging to 
our sample have been taken from White et al. (1997). Fig- 
ure [10] shows the dependence of the radio luminosity of the 
central source, Lr, on the cooling luminosity ^jf, i.e. the X- 
ray luminosity within the cooling volume (we note that, even 
though White et al. 1997 quote their values in an 57o ~ 1, 
ho = 0.5 cosmology, converting these measurements to the 
adopted ACDM cosmology would merely correspond to a 
vertical shift of the ^jf vs Lr relation which does not affect 
our analysis and conclusions). Given the paucity of sources 
and the fact that they did not show any different behaviour, 
in this case the best-fit value for vs Lr is only given 
for the whole sample of clusters which host a central ra- 
dio source, independent of whether this source presents an 
extended or a point-like structure. The plot of the classi- 
cal mass deposition rate versus the radio luminosity shows a 
clear correlation (^ oc ijj with 6 = 0.25±0.02 cfr. Table 1) 
which is in agreement with feedback models (e.g. Churazov 
et al. 2003) for the fueling of cluster radio sources. 



6 CLUSTERS WITH CENTRAL RADIO 
EMISSION: UNDER-LUMINOUS OR 
OVERHEATED? 

In Section 5.2 we have seen how clusters that host a radio 
source endowed with an extended structure show a depar- 
ture from the typical Lx — kT relation found for the cluster 
population as a whole (cfr. Figure [7^), a departure that 



is particularly remarkable in the case of low-mass systems. 
What is not yet clear from our analysis is if extended radio 
sources cause an under-luminosity or an overheating of the 
associated cluster. If it is indeed a temperature increment 
caused by the central radio galaxy, one may expect the cor- 
responding heat excess to be correlated with the radio power 
of the central source. 

To investigate this possibility, we have considered the 
quantity E^cq = 3NkAT/2 = 3M^^skAT/2fimu which is 
the energy required to heat the ICM from the predicted 
to the measured temperature. In the above expression, A'^ 
is the total number of particles, Mgas the gas mass of the 
cluster and AT the temperature increment of clusters which 
host an extended radio source with respect to the Lx — kT 
relation obeyed by clusters with unresolved central sources. 
The gas mass was computed from the total mass of the clus- 
ter by using a constant gas mass fraction of 0.2. We have 
plotted the required energy as a function of radio power in 
Fig. 111! The radio luminosity and the heat input needed 
to produce the observed temperature increment in clusters 
hosting extended radio sources appear to be correlated, al- 
though with a large scatter. This correlation favours a model 
where the temperature increase is caused by radio galaxy- 
induced heating. As noted by Croston et al. (2005) who per- 
formed a similar analysis for groups of galaxies, the large 
scatter in this plot is not surprising given that there are 
many unknown factors such as the age of the source, its 
history and size. 

The issue of temperature increment can be investigated 
further by performing Monte Carlo simulations to be com- 
pared with the X-ray and radio luminosity distributions of 
the present sample. Possible biases due to the joint effects 
of the X-ray and radio selection functions (highlighted by 
the solid curves in the two panels of Figure I12|l , have been 
removed by limiting our analysis to four contiguous redshift 
regions which are complete both in X-ray luminosity and in 
1.4 GHz radio power. In order to maximize the number of 
sources available for statistical analyses, these four regions 
in the z — Lx — Lr space have been identified as follows: 

1) 2 ^ 0.01; Lx ^ 0.0015 ■ 10^"; Lj, ^ lO^^'^ 

2) 0.01 < 2 < 0.05; Lx ^ 0.17 ■ 10^^; Lr ^ 10^^ ^ 

3) 0.05 < 2 < 0.1; Lx > 0.708 ■ 10''''; Lr ^ 10^^'^ 



4) O.K 2 ^ 0.2; Lx > 2.51 ■ 10''''; Lj? > 10^ 



(5) 



- where X-ray luminosities are in erg/s and radio luminosi- 
ties in W/Hz/sr - and are those enclosed within the dashed 
lines in Figure [T^ 91 clusters in our dataset fulfill the above 
requirements (crosses in Figure I12|l , out of which 45 host 
sub-structured/extended radio sources (filled green squares 
with crosses on top in Figure [T^J. 

Simulated samples in the z — Lx — Lr space have been 
obtained by assuming that the cluster X-ray luminosity and 
radio luminosity of the central source are not correlated with 
each other. Each simulated object in the four regions consid- 
ered in our analysis was then given a triplet of values z, Lx 
and Lr, where the X-ray luminosity was assigned according 
to the Bohringer at al. (2002) cluster luminosity function as 
derived from the REFLEX cluster survey once their sample 
was corrected both for missing flux and flux error: 



$(Lx)dix = $0 (3^) '^xp (--j^ 



dLx 



(6) 
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Figure 12. Distribution of X-ray (left-liand panel) and radio (right-hand panel) luminosities as a function of redshift for the 145 
NORAS-I-REFLEX clusters endowed with a central radio-active component. The solid lines show the minimum luminosities probed by 
the X-ray and 1.4 GHz surveys corresponding respectively to the limiting fluxes of 3 ■ 10^'^ erg s"^ cm"'^ and 3 mjy. Open (red) circles 
correspond to point-like radio sources, while filled (green) squares identify sub-structured objects. Crosses highlight those systems (91, 
out of which 45 show signatures of of extended emission) which fulfill the luminosity requirements JSJ adopted for the Monte Carlo 
simulations (dashed horizontal lines both in the left-hand and right-hand panels) and that therefore have been included in the analysis 
presented in this section (sec text for details). 
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Figure 13. Fractional luminosity distributions of X-ray selected clusters inhabited by a central radio-active AGN. The solid points with 
associated errorbars correspond to the 91 sources extracted from the sample presented in this work by requiring the conditions l[5]l to be 
fulfilled. Solid lines represent the results of the Monte Carlo simulations. 



with a = 1.69, $o = 1-07 ■ 10~^/ig,5o Mpc"^ and 
L» — 8.36 ■ 10''''/Iq 5Q erg/s, where /lo.so ~ ho/Q.b. 



proved to provide a very good fit to the local, AGN-fuelled, 
radio population (see e.g. Magliocchetti et al. 2002): 



Radio luminosities were instead drawn from the Dun- 
lop & Peacock (1990) radio luminosity function for steep- 
spectrum sources (model 7 for pure luminosity evolution). 



Lc{z) 



(7) 
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Figure 14. The same as in Figure [T3l but for the sub-population of extended radio sources (45 objects). 



where a = 0.69, /3 = 2.17, *o = lO"'''''' Mpc"^ 
(Alogj^Ql/ji)^^ and where the evolving break luminosity 
Lc{z) can be expressed as logiQLc{z) = 26.22+1.26z — 0.26z^ 
(in Hz~^ units). Since the above expression was derived 
for an Slo = 1, /lo = 0.5 (EdS) universe, following Dunlop & 
Peacock (1990) we have rewritten equation (O in the con- 
cordance ACDM cosmology by using the relation: 

*l(ifl,l,^)^ =*2(Lfl,2,^)^, (8) 

where 1 and 2 refer to the old and new cosmology and dV/dz 
is the volume element. 

We performed 1000 Monte Carlo runs, where each run was 
set so to give the same number of sources per redshift inter- 
val with luminosities above the limits expressed in ((5} as in 
the real dataset. 

The luminosity distributions of X-ray-selected clusters 
inhabited by a radio- loud AGN are shown in Figures [T51 and 
1141 which provide the fractional number of clusters per X- 
ray (left-hand panels) and 1.4 GHz luminosity (right-hand 
panels) for both real data (points with associated errorbars) 
and simulated sample (solid lines. Note that, due to the high 
number of realizations, these values are virtually error-free) . 

Two important conclusions can be drawn from the dis- 
tributions shown in Figures [T3l and [141 . First, the X-ray lu- 
minosity distribution of clusters with a central radio-active 
counterpart is virtually identical to that featured by the 
whole population of X-ray selected clusters brighter than 
the same flux limit. In other words, the presence of a cen- 
tral radio- loud component, however powerful, does not affect 
the cluster X-ray luminosity. This result is found regardless 
to whether the radio AGN shows an extended or point-like 
structure. If we then combine the above flnding with what 
emerged from Sections 5.2 and 5.3, we can conclude that - 
mostly in the case of small mass systems - the interaction 
of extended radio sources with the surrounding environment 
causes a remarkable overheating of the available gas, heat- 
ing which produces the steepening of the Lx — kT relation 



observed in Figure [7^. The contribution of AGN heating 
from a central radio source is less important in the case 
of point-like sources which are not observed to permeate the 
ICM and ii) - independent of the type of radio object - in 
the case of high mass systems which require additional forms 
of energy to balance radiative loss. 

This is in agreement with the results displayed in Fig- 
ure [11] Physically, it means that the radio sources displace 
only a small fraction of the X-ray luminous ICM so that 
there is no noticeable luminosity deficit in clusters with ra- 
dio structure. However, simulations show that if a signifi- 
cant fraction of the jet energy is transferred into potential 
energy thus leading to an expansion of the cluster atmo- 
sphere, also the luminosity of the ICM decreases (see Heinz 
et al. 2006) . Such simulations also suggest that a significant 
fraction of the jet energy is thermalised in the ICM on rel- 
atively short time scales. On long time scales, however, the 
energy injected into the ICM by radio galaxies must end up 
mostly as potential energy by the virial theorem, thus im- 
plying long-term temperature increments to be small. This 
indicates that our analysis is unveiling the short-term effects 
of radio-loud AGN heating. The extended radio structures 
in fact originate from recent AGN activity and are found in 
correspondence with ongoing heating in the ICM. On the 
other hand, those clusters with no extended structures have 
not had a recent (i.e. longer than the life time of a radio lobe) 
outbreak of radio-loud AGN activity and therefore show no 
signs of recent temperature increments. 

While the X-ray luminosity of a typical cluster remains 
unaffected by the presence of a central radio-emitting com- 
ponent, our simulations clearly show that this is not the 
fate for the 1.4 GHz luminosity of a radio galaxy set at 
the bottom of the potential well of a cluster. In fact, the 
right-hand panels of Figures [13] and [14] highlight that there 
is a remarkable deficit of low- luminosity (10^^ W/Hz/sr < 
Ljf, < 10^^ W/Hz/sr) radio sources when compared to the 
whole Fr 3 mjy population. At the same time, radio-loud 
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AGN within the cluster centres seem to favour luminosities 
^ 10^** W/Hz/sr (~ 3cr discrepancy from the results for the 
whole 1.4 GHz radio population). This different behaviour 
between the radio- luminosity distribution of real sources and 
of simulated data sets is found regardless of whether the cen- 
tral radio source presents an extended structure or not. This 
result implies a difference between the luminosity function 
of radio sources sitting in the cluster centres and that of 
the whole radio population, in the sense that the former 
one is much flatter at all luminosities L_r < 10^^ W/Hz/sr. 
Indeed, Figure 4 of Best et al. (2007) clearly shows the dif- 
ferent behaviour of the luminosity functions as derived for 
their sample of BCGs and in the more general case of radio 
sources associated to SDSS galaxies of different mass, even 
though the authors dismiss the evidence as 'tentative'. 

The fact that low- luminosity radio sources are under- 
represented in clusters is most likely due to the 'strangling' 
effect caused by the ovcrdonse central regions of massive 
systems on radio objects that are not powerful enough to 
expand throughout the surrounding medium. On the other 
hand, if we combine the results presented in this section 
with the analysis performed in Section 5.1, we can con- 
clude that extended structures and high radio luminosities 
are more likely to occur in sources that are located in rich 
environments such as groups or clusters of galaxies. This is 
likely caused by a tight interplay between the intracluster 
medium and the energy released by a central radio source, 
which boosts the radio luminosity of the latter one. Indeed, 
in agreement with our results, Barthel & Arnaud (1996) find 
that for a fixed jet kinetic power, radio luminosities can be 
higher in clusters due to the confining effect of the dense 
ICM which reduces losses due to adiabatic expansion in the 
radio lobes, therefore enhancing the radio synchrotron emis- 
sion. 



7 CONCLUSIONS 

Wo have combined the REFLEX and NORAS cluster sur- 
veys with the NVSS dataset to provide a list of X-ray se- 
lected clusters brighter than 3 • 10~^^ erg s~^cm~^ which 
host in their centres (dist < 1.5%rvir) a radio source brighter 
than 3 mjy. Out of 550 clusters set in the whole sky north 
of 6 = —40°, 148 systems (corresponding to ~ 27 per cent 
of the total sample) show signatures of central radio emis- 
sion. We find that three systems were imaged both in RE- 
FLEX and NORAS in the overlapping region 0° < 5 < 2.5° 
between the two surveys. By then removing double identi- 
fications we end up with a sample of 145 clusters that host 
a central radio source with i<i.4GHz ^ 3 mJy. Visual inves- 
tigations of radio maps show that 61 of the radio objects 
(i.e. ^ 42 per cent of the sample) associated with cluster 
centres possess extended radio emission in the form of elon- 
gated blobs, double or triple/multiple structures. Archival 
data (mainly taken from the BAX and VIZIER databases) 
have provided extra information such as temperatures and 
velocity dispersions for some of the systems in our sample. 

The main conclusions of our work can be summarized 
as follows: 

(1) Although the fraction of clusters that host a cen- 
tral radio source brighter than Lr = 10^^ W/Hz/sr is ap- 
proximately constant with cluster mass and equal to 20% 



(in agreement with the results of Best et al. 2007), we 
find that 11 out of 12 (i.e. 92 per cent) local/low-mass 
systems present central radio emission with luminosities 
Lr ^ lO'^'^ W/Hz/sr. This suggests that, either, radio 
sources preferentially inhabit low-mass clusters, or - more 
likely - that by going deep enough in radio flux, most of 
X-ray selected clusters will be found to host a central radio- 
loud AGN. 

(2) The luminosity-temperature relation of clusters that 
host a central radio source follows that found for the whole 
cluster population {Lx oc T^'^ - e.g. White, Jones & Forman 
1997; Arnaud & Evrard 1999; Popesso et al. 2005). However, 
there is a significant discrepancy in the Lx vs kT relation 
between unresolved and extended radio sources. The latter 
population shows the much steeper trend Lx oc T^, while 
for unresolved structures one finds Lx oc T^'^, which is close 
to the self-similar result. The difference between the Lx vs 
kT trend for point-like and extended sources becomes par- 
ticularly noticeable for temperatures kT < 3 KeV, where 
almost all systems with extended radio emission lie below 
those associated with unresolved radio structures. 

(3) The radio luminosities of the central sources show a 
steep {Lr (xT^) dependence on the cluster temperature for 
kT < 3 KeV. Such a correlation is lost in more massive 
systems. 

(4) The X-ray luminosity vs velocity dispersion and the 
temperature vs a relations of our sample are in agreement 
with those found in the literature (Lx oc a"' *; kT oc a^; e.g. 
Popesso et al. 2005; Madhavi & Gheller 2001; Osmond & 
Ponman 2004) . No significant differences are found between 
point-like and extended sources, except for a mild preference 
for the latter population to appear in higher velocity sys- 
tems. We remark though that the available cr measurements 
do not allow us to probe the very low-mass regime where the 
different behaviour between point-like and extended sources 
could be more prominent. 

(5) The mechanical luminosity provided by the central 
radio source is found to balance radiative losses only in small 
(i.e. a < 400 km/s) systems. 

(6) It seems plausible that heating by low-power (FR I) 
radio galaxies can explain the absence of cooling flows in 
clusters of galaxies. This is in agreement with conclusions 
by various other authors (e.g. Fabian et al. 2003; Croston 
et al. 2005), who concentrated their analysis on groups of 
galaxies. However, for the largest clusters, the energy re- 
leased by the radio sources may not be sufficient to balance 
radiative losses. 

(7) Monte Carlo simulations show that the X-ray lumi- 
nosity of a cluster is not affected by the presence of a central 
radio source, however powerful. The luminosity distribution 
of the present sample of radio-identifled systems is the same 
as that of the whole cluster population. Combining this re- 
sult with what we found in point (2), we can conclude that 
the presence of radio sources with an extended structure is 
responsible for the over-heating of the intracluster gas. The 
importance of such a heating dramatically increases in low- 
mass systems. 

(8) The luminosity distribution of radio galaxies sitting 
in cluster centres is very different from that of the total 
radio population. In fact, in the former case we flnd that 
low luminosities {Lr ~ 10^^ — 10^^ W/Hz/sr) are depressed, 
while higher luminosities are strongly boosted. The net effect 
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on the radio luminosity function is a flattening in the whole 
luminosity range Lr < 10^* W/Hz/sr. 

Our results support a strong interaction between AGN 
radio emission and the ICM. Radio sources are present in 
the centres of almost all local/small meiss systems and we 
expect the majority of clusters to exhibit central radio emis- 
sion with powers Lr > 10^° W/Hz/sr. The ICM strongly 
affects the luminosity of the radio source that sits in the 
cluster centre: low-luminosity objects are 'strangled' by the 
overdense central medium, while bright ones have their lu- 
minosities boosted by the interaction with the surrounding 
gas (see e.g. Barthel & Arnaud 1996). Moreover, the exis- 
tence of a central radio source - especially if it exhibits an 
extended structure and resides in small-mass systems - leads 
to a significant heating of the ICM. 

Furthermore, radio-loud AGN are observed to provide 
enough mechanical energy to balance radiative losses. The 
effect of AGN heating from a central radio source becomes 
less important t) in the case of point-like sources which are 
not observed to permeate the ICM and ti) - independent of 
the type of radio object - for high mass systems. Clusters 
hotter than kT < 3 KeV, do not obey the tight Lr vs kT re- 
lation observed for small systems, the effect on the thermal 
state of the cluster as provoked by the presence of a radio- 
active AGN is independent of the capabilities of such an 
object to permeate the surrounding medium (no difference 
in the Lx — kT relation between extended and point-like 
sources) and AGN heating is found to be insufficient to bal- 
ance radiative losses in the cluster. In this case, a different 
source of heating such as thermal conduction (e.g. Narayan 
& Medvev 2001; Voigt & Fabian 2004) has to be invoked. 

The data suggest a smooth transition between the 
radio-AGN heating mode and the thermal conduction 
mode. Indeed, double heating models - first developed by 
Ruszkowski & Begelman (2002) - have been proved to pro- 
vide a good agreement with the observed cluster properties 
(see e.g. Briiggen et al. 2003; Hoeft & Briiggen 2004; Roy- 
chowdury et al. 2005; Fujita & Suzuki 2005). Within this 
framework, Hoeft & Briiggen (2004) have shown the rela- 
tive importance of heating due to thermal conduction to 
increase with cluster mass, in excellent agreement with our 
findings. 

Obviously, it would be desirable to push this kind of 
analysis to fainter fiuxes (both radio and X-ray) and higher 
redshifts to investigate whether some of our findings only 
hold in the local universe or are more general properties of 
the cluster population. Forthcoming surveys like COSMOS 
(see e.g. Finoguenov et al. 2007) are expected to provide 
such answers. 
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